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ABSTRACT 

Context. Gamma-ray burst (GRB) emission is believed to originate in highly relativistic fireballs. 

Aims. Currently, only lower limits were securely set to the initial fireball Lorentz factor Tq. We aim to provide a direct measure of Tq. 
Methods. The early-time afterglow light curve carries information about To, which determines the time of the afterglow peak. We 
have obtained early observations of the near-infrared afterglows of GRB 060418 and GRB 060607A with the REM robotic telescope. 
Results. For both events, the afterglow peak could be clearly singled out, allowing a firm determination of the fireball Lorentz of 
r ~ 400, fully confirming the highly relativistic nature of GRB fireballs. The deceleration radius was inferred to be W dcc ss 10 17 cm. 
This is much larger than the internal shocks radius (believed to power the prompt emission), thus providing further evidence for a 
different origin of the prompt and afterglow stages of the GRB. 

Key words. Gamma rays: bursts - Relativity 



1. Introduction 

The early stages of gamma-ray burst (GRB) afterglow light 
curves display a rich variety of phenomena at all wave- 
lengths and contain significant information which may allow 
determining the physical properti es of the emitting f ireball. 
The launch of the Swift satellite dGehrels et alj l2004l) . com- 
bined with the development of fast-slewing ground-based tele- 
scopes, has hugely improved the sampling of early GRB af- 
terglow light curves. To date there are several published early 
optical/near-infrared (NIR) afterglow light curves that cannot 
be fitted with a simple power-law, testifying that the detec- 
tions started before the afterglow began its regular decay. To 
explain the early b ehaviour of m ost of them (GRB 990123 : 
lAkerlofetalJ 119991: GR B 0210 04: iKobavashi & Zhang] 120031: 
GRB 021211: |Fbx_eta]J llQQl; |Li_eta]J |2003j GRB041219A: 
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Fan et al.l206HlVestrand_et al.l2005t GRB 50730:rPandev etal.1 
20061 GRB 05 0801: iRoming et all | 2006t iRvkoff et al l 120061: 
GRB 50820A: IVestrand e t al.1 120061: G RB 050904: iBoer et al I 
12001 IWeietall I2006t GRB 06020 6: [Wozniak et al.1 12 006: 
iMonfardini eTail l2006t GRB 060210: IStanek et all 120071) dif- 
ferent mechanisms have been proposed, i.e. reverberation of 
the prompt emission radiation, reverse and refreshed shocks 
and/or energy injection a nd, lately, large angle emission 
dPanaitescu & Kumdl2007l) . 

Since many processes work in the early afterglow, it is often 
difficult to model them well enough to be able to determine the 
fireball characteristics. The simplest case is a light curve shaped 
by the forward shoc k only, as could be the case for GRB 03 0418 
(IRvkoff et al.ll2004l) and GRB 060124 dRomano et alj|2006l) . but 
the lack of a measured redshift in the former case and the poor 
sampling in the latter prevented to derive firm conclusions. This 
case is particularly interesting because, while the late-time light 
curve is independent of the initial conditions (the self-similar 
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Fig. 1. NIR and X-ray light curves of GRB 060418. The dotted 
lines show the models of the NIR data using the smoothly broken 
power law (see Sect. 13.11 ). while the dashed line shows the best- 
fit to the X-ray data. The vertical lines mark the epochs at which 
we computed the SED (Fig. [3}. 

solution), the time at which the afterglow peaks depends on the 
original fireball Lorentz factor T, thus all owing a direct mea - 
surement of this fundamental parameter (ISari & Piranl [1 999). 
The short variability timescales, coupled with the nonthermal 
GRB spectra, indeed imply th at the sources emi t ting GRBs have 



a highly relativistic motion (Rudermanl 19751 [Fenimor e et al.l 
1 1993] iPiranl 120001: iLithwick & Saril 1200 ll) . to avoid suppres- 
sion of the high-energy photons due to pair production. This 
argument, however, can only set a lower limit to the fireball 
Lorentz factor. Late-time measure ments (weeks to months after 
the G RB) have shown T ~ a few dFrail et al-ll 19971 iTavlor et al. I 
120051) . but a direct measure of the initial value (when Y is ex- 
pected to be ~ 100 or more) is still lacking. 

We present here the NIR early light curves of the 
GRB 060418 and GRB 0606 07 A afterglows observed with the 
REM robotic telescoped dZerbi etalj l2001t IChincarini et all 
120031) located in La Silla (Chile). These afterglows show the on- 
set of the afterglow and its decay at NIR wavelengths as simply 
predicted by the fireball forward shock model, without the pres- 
ence of flares or other peculiar features. 

2. Data 

GRB 060418 and GRB 060607 A were detected by Swift 
at 03:06:08 UT (iFalcone et all l2006ah and 05:12:13UT 
dZiaeepour et al.l 120061) . respectively. The BAT light curve 
of the former (Tgp - 52 + 1 s) showed three overlapping 
peaks (Cummings et al. 2006). For the latter, the light curve 
is dominated b y a double-peaked structure with a duration 
T 90 = 100 ± 5 s dTueller et al.ll2006l) . The Swift XRT started ob- 
serving the fields 78 and 65 s after the trigger, respectively. For 
both bursts a bright, previously uncatalogued fading source was 
singled out. The XRT light curves show prominent flaring activ- 
ity, superposed to the regular afterglow decay. For GRB 060418 
a prominent peak, also visible in the BAT data, wa s observed 
at about 128 s after the trigger dFalcone et alj|200 6b) while the 



Fig. 2. H-band and X-ray light curves of GRB 060607 A. The 
dotted line shows the model to the NIR data using the smoothly 
broken power law (see Sect. 13. U . 

Table 1. Best fit values of the light curves of the first 
hour of observations for GRB 060418 and the first 1000 s for 
GRB 060607 A (lcr errors), using the smoothly broken power- 
law reported in Sect. O with f pea k = fb(-ffi-/«d) 1/Wffd ~" r)] . The 
relatively large x 1 of the fit results from small-scale irregularities 
present throughout the light curve (see Figs. [T]and|2]l. 
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XRT light c urve of GRB 0606 07 A was characterised by at least 
three flares dPage et aT]|2006l) . UVOT promptly detected bright 
optical counterparts for both events. Th e redshift is z — 1 -489 
for G RB 060418 (iDupree et al.l 120061: IVreeswiik & Jaunsen I 
120061) and z = 3.082 for GRB 060607 A dLedoux et alJl2006h . 

The REM telescope reacted promptly to both GCN alerts 
and began observing the field of GRB 060418 64 s after the 
burst (39 s after the reception of the alert) and the field of 
GRB 060607 A 59 s after the burst (41 s after the reception of 
the alert). REM hosts REMIR, an infrared imaging camera that 
operates in the range 1.0-2.3 /im (z'JHK'), and ROSS, an op- 
tical imager and slitless spec trograph. For both targets a bright 
NIR source was identified dCovino et al.ll2006allbl) . ROSS could 
not observe these two GRBs due to maintenance work. REM fol- 
lowed the two events down to the sensitivity limits adopting two 
different observing strategies. In the former case multi-colour 
z' JHK-band observations were carried out while in the latter 
case only the H filter was used in order to get a denser sampling 
of the light curve in just one band. All our photometric data are 
reported in Tables [2] and [3] 

3. Results and Discussion 

3. 1 . Light curve modelling 

Figures Q] (GRB 0604 1 8) and [2] (GRB 060607 A) show the NIR 
and X-ray light curves of the two afterglows. The X-ray data 
have been taken with the Swift XRT. The analysis was carried 
out adopting the standard pipeline. Pile-up corrections both in 
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WT and PC mode were taken into account when required. The 
NIR light curves of the two events show a remarkable similarity. 
Both present an initial sharp rise, peaking at 100-200 s after 
the burst. The NIR flux of GRB 060418 decays afterwards as a 
regular power law. The NIR light curve of GRB 060607A shows 
a similar, smooth behaviour up to ~ 1000 s after the trigger, 
followed by a rebrightening lasting ~ 2000 s. 

To quantitatively evaluate the peak time, we fitted 
the NIR light curves using a smoothly broken power-law 
dBeuermann et all 1 19991) : F(t) = F /[(t/t b ) Ka ' + (t/t b Y a "] 1/K , 
where Fq is a normalisation constant, a T (d) is the slope of the 
rise (decay) phase (a, < 0) and a: is a smoothness param- 
eter. The time at which the curve reaches its maximum is 
fpeak = f b (-Q'r/Q'd) 1/Wl,d ^ )] - We obtain for GRB 060418 and 
GRB060607A peak times of 153 + 10 and 180 + 6 s, respec- 
tively. Any other suitable functional forms provides comparable 
results. The complete set of fit results is reported in Table [TJ 

As for many other GRBs observed by Swift, the early X-ray 
light curves of both events show several, intense flares superim- 
posed on the power-law decay (IChincarini et al 1 I2007T) . In par- 
ticular, for GRB 060418 a bright flare was active between ~ 115 
and 185 s. Excluding flaring times, the decay is then described 
by a power law with decay slope ax = 1.42 + 0.03. The X- 
ray light curve of GRB 060607A is more complex and presents 
two large flares within the first 400 s. After that the flux den- 
sity decreases following a shallow power law (with small-scale 
variability), until steepening sharply at t ~ 10 4 s. 

By comparing the X-ray and NIR light curves of both bursts, 
it is apparent that the flaring activity, if any, is much weaker at 
NIR frequencies. It is thus likely that the afterglow peak, vis- 
ible in the NIR, is hidden in the X-ray region. In the case of 
GRB 060418, thanks to our multicolour data, this can be con- 
firmed by considering the NIR/X-ray spectral energy distribu- 
tion (SED; Fig. g). At f = 800 s, when both the NIR and X- 
ray light curves are decaying regularly, the SED is described 
by a synchrotron spectrum, with the cooling frequency lying at 
Vcooi ~ 3 x 10 15 Hz. The spectral slopes at NIR and X-ray wave- 
lengths are /? NIR = 0.65 +0.06 and /3 X = 1.03 +0.04, respectively. 
A small amount of extinction in the host (Ay = 0.10 mag as- 
suming the SMC extinction curve) was considered to make /?nir 
consistent with the synchrotron expectation ^Snir = fix -0.5. The 
situation is clearly different at early times, when X-ray flaring is 
active: the X-ray component is softer (fix = 1.27 + 0.03) and 
much brighter than the extrapolation of the NIR emission (even 
including the dust correction). This implies a different origin for 
the X-ray emission. For GRB 060607 A, the light curve is com- 
plicated by several flares. As we lack multicolour data, a detailed 
analysis is not possible. 



3.2. Determination of the Lorentz factor F 

Our spectral and temporal analysis agrees with the interpretation 
of the NIR afterglow light curves as corresponding to the af- 
terglow onset, as predicted by the f ireball forward shock model 
(Sari &Piranlll999l [Meszaros 20061) . According to the fits, the 
light curves of the two afterglows peak at a time f pea k > Tgo 
as expected in the impulsive regime outflow ('thin shell' case). 
In this scenario the quantity f pea k/(l + z) corresponds to the de- 
celeration timescale fd ec ~ ^dec/(2cF^ ec ), where /?d ec is the de- 
celeration radius, c is the speed of light and Fd ec is the fireball 
Lorentz factor at fd ec- It is therefore possible to estimate Tdec 
dSari & Piranl 19991). which is expected to be half of th e ini- 
tial value r dPanaitesc u"& Kumarll2000l: lMeszarosll2006l) . For 
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Fig. 3. SED of the GRB 060418 afterglow at 130 and 800 s after 
the trigger. NIR data have been corrected for Galactic and host 
extinction (Ay — 0.74 and 0.10 mag, respectively). 

a homogeneous surrounding medium with particle density n, we 
have 



F(fpeak) 
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Jy — iu ii r ,53 erg is the isotropic-equivalent en- 
ergy released by the GRB in gamma rays, n = «o cm -3 , 
fpeak,2 = fpeak/(100 s), t] = 0.2 770.2 is the radiative efficiency 
and OT p is the proton mass. We use E y — 9 x 10 52 erg for 
GRB 060418 dGolenetskii et al.ll2006l) a nd E y ~ 1.1 x 10 53 erg 
for GRB 060607A dTueller et al.1 l2006h . Substituting the mea- 
sured quantitie s and normalising to the typical values n — 1 cirT 3 
and i] = 0.2 (Bloo m et al]|2003l) . we infer for both bursts Tq ~ 
400 (?7o.2io) _1 ■ This value is very weakly dependent on the un- 
known parameters n and rj, and therefore provides a robust de- 
termination of Tq. 

One possible caveat concerns the density profile of the 
surrou nding medium, which is still uncertain. Chevalie r & Lil 
(|2000|) proposed a wind-shaped profile n{r) = Ar~ 2 , where A 
is a constant, as expected around a massive star. In principle, 
the temporal and spectral properties of the afterglow can be used 
to distinguish among the two cases. For v < v coo i, before fdec, 
the flux is expected to rise as f 3 for the homogeneous (ISM) 
case, while it evolves slower than t 1 ^ in a wind environment. 
The measured values for GRB 060418 and GR B 060607A seem 
therefore consistent with the ISM case (see also lJin &~F an 2007). 
For GRB 060418, however, no closure relation is satisfied for 
t > fpeak- In fact, irrespective of the density profile, the tem- 
poral and spectral slopes in the X-ray region (which lies above 
Vcooi) should verify ax = 3/3x/2 - 1/2, which is clearly not 
the case given the observed values. Such discrepancy has been 
observed in many other examples and may be due to radiative 
losses, o r varying equipartition parameters, o r Compton losses 
(see e.g. iPanaitescul 120051: IZhang et alll2006h . A detailed treat- 
ment of these effects is beyond the scope of this paper. We how- 
ever provide an estimate of To in the case of a wind environment: 
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Using the same values as above for E and z, and assuming A* - 
A/(3 x 10 35 cm- 1 ) = 1, we find T ~ 150 for both bursts. 

The determination of f pea k is in principle affected by the 
choice of the time origin tp, which we set to the B AT trigger time 
dLazzati & Begelmanll200rl lOuimbv et al.ll2006l) . These authors 
have however shown that this effect is small, and mostly affects 
the rise and decay slopes rather than the peak time (especially 
when fp ea k is larger than the burst duration, as for GRB 060418 
and GRB 060607 A). The measurement of To is thus not very 
sensitive to the exact choice of to. 

For both bursts, we could not detect any reverse shock emis- 
sion. The lack of such flashes has already been noticed previ- 
ously in a set of Sw ift bursts with prompt UVOT observations 
dRoming et al.l 120061) . Among the many possible mechanisms 
to explain the lack of this component, strong suppression (or 
even total lack) of reverse shock emission is n aturally expected 
if the outflow is Poyntin g-flux dominated (Fan et al. 2004; 
IZhang & Kobayashil 120051) . Nevertheless, iJin & Fad ( 120071) 
showed that for GRB 0604 18 and GRB 060607 A the reverse 
shock emission might be too weak to be detected. 

Finally, the optical afterglow of GRB 030418 showed a 
pe aked light cur v e simi lar to our cases, and has been explained 
by iRvkoffet al.l d2004) as due to decreasing extinction. Our 
peaks are too sharp to match this interpretation. Furthermore, 
a very large extinction would be implied, which contrasts with 
the UVOT detect ion in the rest-frame UV dSchady et al.ll2006t 
IOatesetal.ll2006l) . The peak of GRB 030418 might have been 
the afterglow onset. Setting z — 1, we get Fo ~ 100 using Eq. ([T). 

4. Conclusions 

The REM discovery of the afterglow onset has demonstrated 
once again the richness and variety of physical processes oc- 
curring in the early afterglow stages. The very fast response 
observations presented here provide crucial information on the 
GRB fireball parameters, most importantly its initial Lorentz 
factor. This is the first time that F(f pea k) is directly measured 
from the observations of a GRB. The measured To value is 
well within the range (50 ^ rp ^ 1 000) envisaged b^ 



the standard fireball model (Pira 



ISoderberg & Ramirez-Ruizll2002t Meszaros 200 61). It is also in 
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agreement with existing measured lower limits ( Lithwick & Saril 
120011: IZhang et al.l l2006). Albeit the values for Fo are compara- 
ble for GRB 060418 and GRB 060607 A, as noted in the intro- 
duction prompt optical detections have shown very different be- 
haviours, so this value should not be taken as common to every 
GRB. 

Using Tp = 400 we can also derive other fundamental 
quantities characterising the fireball of the two bursts. In par- 
ticular, the isotropic-equivalent baryonic load of the fireball is 
Mfb = £/(rpc 2 ) 7 x 1O~ 4 M , and the deceleration radius is 
R dec « 2cfp eak [T(fpeak)] 2 /(l + z) ~ 10 17 cm. This is much larger 
than the scale of ~ 10 15 cm where th e internal shocks are be- 
lieved to power the pro mpt emission (M eszaros &Reeslll997t 
Rees & Meszaros 1994), thus providing further evidence for a 
different origin of the prompt and afterglow stages of the GRB. 
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Table 2. Observation log for GRB 060418. The time t indicates 
t he BAT trigger tim e 2006 April 18.12926 UT. 
* lNvsewander et al.l d2006l) . 



Mean time 


t-t 


Exp. time 


Filter 


Magnitude 


(UT) 


(s) 


(s) 








Apr 18.12974 


40 


5 


z 


15.3: 


±0.3* 


Apr 18.13371 


479 


100 


zf 


13.26 


±0.06 


Apr 18.14215 


1114 


100 


t! 


14.70 


±0.09 


Apr 18.15722 


2416 


100 


z' 


15.58 


±0.12 


Apr 18.13264 


292 


100 


J 


12.33 


±0.05 


Apr 18.13975 


906 


150 


.1 


13.82 


±0.06 


Apr 18.15323 


2071 


300 


J 


14.94 


±0.07 


Apr 18.13014 


76 


10 


H 


11.98 


±0.16 


Apr 18.13032 


92 


10 


H 


11.46 


±0.07 


Apr 18.13050 


107 


10 


H 


11.00 


±0.19 


Apr 18.13067 


122 


10 


H 


10.92 


±0.12 


Apr 18.13085 


137 


10 


H 


11.01 


±0.02 


Apr 18.13538 


529 


150 


H 


12.24 


±0.03 


Apr 18.14527 


1383 


150 


H 


13.57 


±0.04 


Apr 18.16008 


2663 


150 


H 


14.89 


±0.22 


Apr 18.18578 


4883 


150 


H 


15.48 


±0.10 


Apr 18.19649 


5809 


150 


H 


15.40 


±0.12 


Apr 18.20720 


6734 


150 


H 


15.60 


±0.10 


Apr 18.21791 


7659 


150 


H 


15.83 


±0.15 


Apr 18.13121 


168 


10 


K' 


10.23 


±0.04 


Apr 18.13139 


184 


10 


K' 


10.26 


±0.03 


Apr 18.13157 


200 


10 


K' 


10.38 


±0.03 


Apr 18.13175 


215 


10 


K' 


10.48 


±0.05 


Apr 18.13192 


230 


10 


K' 


10.49 


±0.05 


Apr 18.14652 


1491 


300 


K' 


13.14 


±0.09 



Table 3. Observation log for GRB 060607A. The time f indi- 
cates the BAT trigger time, 2006 June 7.21682 UT. 



Mean time 
(UT) 


t-h 
(s) 


Exp. time 

(s) 


Filter 


Magnitude 


Jun 7.21767 


73 


10 


H 


14 


.60 = 


= 0. 


,20 


Jun 7.21794 


97 


10 


H 


13 


.35 i 


= 0. 


.10 


Jun 7.21811 


111 


10 


H 


12 


.91 i 


= 0. 


.07 


Jun 7.21829 


127 


10 


H 


12 


.57 = 


= 0. 


.06 


Jun 7.21848 


143 


10 


H 


12 


.33 = 


= 0. 


.05 


Jun 7.21866 


159 


10 


H 


12 


.01 i 


= 0. 


.05 


Jun 7.21883 


174 


10 


H 


12 


AS-. 


= 0. 


.05 


Jun 7.21902 


190 


10 


H 


12 


.10 = 


= 0. 


.05 


Jun 7.21919 


205 


10 


H 


12 


.07 = 


= 0. 


.05 


Jun 7.21937 


220 


10 


H 


12 


.21 = 


= 0. 


.05 


Jun 7.21955 


236 


10 


H 


12 


.21 = 


= 0. 


.05 


Jun 7.21973 


251 


10 


H 


12 


.31 ±0. 


.05 


Jun 7.21991 


267 


10 


H 


12 


.40 = 


= 0. 


.05 


Jun 7.22008 


282 


10 


H 


12 


.40 = 


= 0. 


.05 


Jun 7.22027 


298 


10 


H 


12 


.55 = 


= 0. 


06 


Jun 7.22045 


314 


10 


H 


12 


.72 = 


= 0. 


.06 


Jun 7.22062 


328 


10 


H 


12 


.59 = 


= 0. 


.06 


Jun 7.22080 


344 


10 


H 


12 


.82 = 


= 0. 


.07 


Jun 7.22098 


359 


10 


H 


12 


.84 = 


fcO. 


.07 


Jun 7.22128 


385 


30 


H 


12 


.95 = 


= 0. 


.05 


Jun 7.22169 


421 


30 


H 


12 


.96 ± 0. 


.05 


Jun 7.22209 


455 


30 


H 


13 


.10 = 


= 0. 


.05 


Jun 7.22251 


492 


30 


H 


13 


.23 = 


= 0. 


.05 


Jun 7.22292 


527 


30 


H 


13 


.25 = 


= 0. 


.05 


Jun 7.22400 


620 


150 


H 


13 


.40 = 


= 0. 


.04 


Jun 7.22606 


798 


150 


H 


13 


.96 = 


= 0. 


06 


Jun 7.22812 


976 


150 


H 


14 


.32 = 


= 0. 


.07 


Jun 7.23033 


1167 


150 


H 


14 


.94 = 


= 0. 


.08 


Jun 7.23223 


1331 


150 


H 


14 


.63 = 


= 0. 


.07 


Jun 7.23427 


1508 


150 


H 


14 


.69 = 


= 0. 


.08 


Jun 7.23635 


1687 


150 


H 


14 


.55 = 


= 0. 


.07 


Jun 7.24047 


2043 


150 


H 


14 


.49 = 


= 0. 


.07 


Jun 7.24253 


2221 


150 


H 


14 


.68 = 


= 0. 


.08 


Jun 7.24459 


2399 


150 


H 


15 


.05 = 


= 0. 


.09 


Jun 7.24767 


2665 


300 


H 


15 


.46 = 


tO. 


.08 


Jun 7.25113 


2964 


300 


H 


15 


.73 = 


= 0. 


.10 


Jun 7.25493 


3293 


300 


H 


15 


.58 = 


= 0. 


.09 


Jun 7.25886 


3632 


150 


H 


15 


.48 = 


= 0. 


.13 


Jun 7.37574 


13730 


6000 


H 


17 


.23 = 


= 0. 


.23 



